Quadrupolar spectra of nuclear spins in strained InGaAs quantum dots 
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Self-assembled quantum dots (QDs) are born out of lattice mismatched ingredients where strain 
plays an indispensable role. Through the electric quadrupolar coupling, strain affects the magnetic 
environment as seen by the nuclear spins. To guide prospective single- QD nuclear magnetic reso- 
nance (NMR) as well as dynamic nuclear spin polarization experiments, an atomistic insight to the 
strain and quadrupolar field distributions is presented. A number of implications of the structural 
and compositional profile of the QD have been identified. A high aspect ratio of the QD geometry 
enhances the quadrupolar interaction. The inclined interfaces introduce biaxiality and the tilting of 
the major quadrupolar principal axis away from the growth axis; the alloy mixing of gallium into 
the QD enhances both of these features while reducing the quadrupolar energy. Regarding the NMR 
spectra, both Faraday and Voigt geometries are investigated, unraveling in the first place the extend 
of inhomogeneous broadening and the appearance of the normally-forbidden transitions. Moreover, 
it is shown that from the main extend of the NMR spectra the alloy mole fraction of a single QD can 
be inferred. By means of the element-resolved NMR intensities it is found that In nuclei has a factor 
of five dominance over those of As. In the presence of an external magnetic field, the borderlines 
between the quadrupolar and Zeeman regimes are extracted as 1.5 T for In and 1.1 T for As nuclei. 
At these values the nuclear spin depolarization rates of the respective nuclei get maximized due to 
the noncollinear secular hyperfine interaction with a resident electron in the QD. 

PACS numbers: 75.75.-c, 76.60.Gv, 76.60.Pc 
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I. INTRODUCTION 

In the coherent control of electron spins within a 
solid state environment, such as in quantum dots (QDs), 
the nuclear spin reservoir acts as the main source of 
decoherenceii"— Recent studies now assure that nuclear 
spin bath can be tamed so as to counter the detrimen- 
tal effect it may have on the carrier spins There- 
fore, in the emerging state of understanding, an alter- 
native is to utilize nuclear spins as a resource;^ for in- 
stance as a natural qubit memory (for an extended re- 
view and references, see, Ref. |10|). In these studies, 
self- assembled QDs have been one of the archetypal sys- 
tems. Within the full gamut of possible combinations 
of self-organized materials;^ the InAs on GaAs system 
stands out due to extensive research efforts devoted over 
more than a decade to their growth and optoelectronic 
characterizations.^^ For the purposes of controlling the 
spin dynamics in InAs QDs, several groups^^"— have 
studied various aspects of the optical orientation^^ of an 
electron spin and its transfer to the nuclear spins via the 
hyperfine interaction leading to the dynamic nuclear spin 
polarization;^ 

One of the remarkable achievements that paved the 
way to the coherent manipulation of nuclear spin dynam- 
ics has been the nuclear magnetic resonance (NMR) of 
a single QD.^^"^^ With further improved precision, po- 
sition selective control of small groups of nuclear spins 
inside a dot is being pursued. This so-called optically 
detected NMR is initiated by the circularly polarized ex- 
citation of a spin-polarized electron spin inside a QD that 
polarizes the nuclear spins through the hyperfine interac- 



tion. The Overhauser field established by the polarized 
nuclei acts back on the electronic system, which can be 
externally measured over an excitonic Zeeman splitting. 
Additionally if an rf magnetic field resonant with nuclear 
spin transitions is incident, it depolarizes some of the 
nuclear spins reducing the Overhauser field, which can 
in turn be detected from a shift in the optical emission 
spectra. Along this line, a recent demonstration utilized 
a sequence of two phase-locked rf pulses to induce coher- 
ent rotations of a targeted group of nuclear spins optically 
pumped to a high polarization degree;^ The advantage 
of this technique is that it enables full coherent control 
over the Bloch sphere and yet on the order of microsec- 
ond time scales. 

These NMR experiments were performed on GaAs 
interface-fiuctuation QDs with the deliberate aim of 
avoiding any strain to keep the resonances narrow 
On the other hand strain is an integral part of self- 
assembled InAs QDs.^^ The existing negative sentiments 
for the strain in the context of electron and nuclear spin 
dynamics have been dramatically reversed by the work 
of Dzhioev and KorenevjS?, Actually, several decades ago 
it was experimentally shown that anisotropic strain in a 
III-V crystal lattice causes local electric field gradients 
(EFG) with which a spin-J nucleus with I > 2 interacts 
because of its quadrupolar moment<^"— This quadrupo- 
lar interaction (QI) splits the nuclear spin degeneracy 
even in the absence of an external magnetic fieldj^i^l 
hence it energetically suppresses the nuclear spin flip 
events, stabilizing the electron spin orientation confined 
in the QD.^^ With this paradigm shift, strain is no longer 
just a nuisance but something to be exploited as a new 
degree of freedom to tailor the magnetic environment of 
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InAs QDs. 

The aim of this work is to offer an atomistic under- 
standing of the interesting physics arising from the co- 
existence of QI together with the dc and rf magnetic 
fields. This is in accord with the current progress of the 
experimental techniques having the goal to address and 
control relatively few number of nuclei within a QDi^i^ 
Starting from the behavior of the spatial variation of the 
strain tensor, we trace the factors that affect QI, and 
identify primarily the high aspect ratio of the QD as a 
factor for enhancing it. Considering different magnetic 
field orientations and the effect of random alloy mix- 
ing i.e., Ina;Gai_3;As QDs, critical insight is gained on 
the vulnerability of each nuclear spin transition in the 
NMR spectrum under inhomogeneous broadening. An 
important complication is that the mole fraction of the 
constituents can vary appreciably from dot to dot within 
the same samplej^^i^^ rendering the in situ compositional 
identification of a targeted QD so far not practical. We 
demonstrate that this can be readily extracted from the 
main span of the NMR spectra. Moreover, our element- 
resolved spectra present crucial information for the la- 
beling of the features in the rather complicated overall 
spectra. For the typical QD parameters reported in rel- 
evant experiments, we extract an effective quadrupolar 
field, Bq of 1.5 T (1.1 T) for the In (As) nuclei, as a bor- 
derline below which QI dominates. The atomistic picture 
that we acquire also provides us the distribution of the 
quadrupolar principal axes within the QD. We make use 
of this information in working out the nuclear spin de- 
polarization due to noncollinear secular hyperfine inter- 
action arising from the tilting of the major quadrupolar 
axis from the optical axis. The depolarization time drops 
to a minimum on the order of an hour at the magnetic 
field coinciding with Bq , as a manifestation of the strong 
competition between the QI and the Zeeman field at this 
value. 

Regarding the organization of the paper: in Sec. II 
we give the theoretical details for the strain, QI, and 
the nuclear spin depolarization; in Sec. Ill some informa- 
tion about our QD structures is given, followed by the 
results grouped into strain and quadrupolar splittings, 
NMR spectra, and the noncollinear secular hyperfine in- 
teraction; in Sec. IV we summarize our conclusions, and 
append a section on the matrix elements for obtaining 
the energy spectra and the rf-initiated nuclear spin tran- 
sitions under an arbitrary EFG. 



II. THEORY 



Although the theoretical ingredients employed in this 
work are not new, an account of the procedures and the 
mathematical models will be helpful for clarifying how 
the results are obtained. 



A. Strain 

From a computational point of view, to identify the 
strain profile, the ionic relaxation of the QD and the host 
matrix atoms is needed. The technique used for this pur- 
pose is molecular statics as implemented in the LAMMPS 
code.'*° Here, the main input is the semi-classical force 
field that governs the atomic interactions. The preferred 
choice for group IV and III-V semiconductors are the 
Abell-Tersoff potentials. We make use of the recent 
parametrization by Powell et al. who fitted their force 
field expression to a large set of cohesive and elastic prop- 
erties obtained from density functional theoryjiSi With 
these tools and for a chosen compositional profile, the re- 
laxed atomic positions become readily available ."^2, Next, 
one needs to extract the strain state, again in an atom- 
istic level as we aim for the EFGs at each nuclear site. 
Among the several possible strain measures, we adopt 
the one proposed by Pryor et al because of its inherent 
compatibility with tetrahedral bonds, as in the present 
case4^ In this measure, the critical task is to correctly 
align an unstrained reference tetrahedron for each lo- 
cal bonding topology. This local strain attains a very 
rapid variation especially for the random alloy mixing of 
In^Gai-ajAs. Therefore, like Bester et al.,^ we average 
the strain tensor; in our case over the twelve neighboring 
same-ionicity sites, which we also apply to the pure InAs 
QDs as well. 



B. Quadrupolar interaction 

In the linear elastic limit, we can express the EFG ten- 
sor components Vij in any orthogonal coordinate framed 
using the computed local strain tensor as: 



dxidxj 



k.l=l 



kl^kU 



(1) 



where S is the fourth-rank gradient elastic tensor4ii^ 
Transforming this tensor expression to Voigt notation in 
the cubic crystallographic frame we get 



(2) 



1^=1 



As the trace of the EFG is unobservablea* it is conve- 
niently set to zero, J^i^a ~^ 0, which leads to 6*11 — 
—25*12. This results in the following explicit relations 
for the EFG tensor components in the mixed Voigt and 
tensor notation 



14, 



Sii 



1 



(3) 

^ xy — ^ yx — 25*4462;^ , (4) 

with the remaining components being obtained by their 
cyclic permutations.^^ The EFG tensor couples to the nu- 
clear quadrupole moment tensor (operator) Q^p through 
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the Hamiltonian^l 
eQ 



6/(2/ 



where I is the dimensionless spin operator, e is the elec- 
tronic charge, / is 9/2 for In, and 3/2 for As and Ga 
nuclei, and Q is the electric quadrupole moment of the 
nucleus. This expression gets simplified in the frame of 
EFG principal axes (V/j = 0, for / ^ J) as 



nc 



e'qQ 
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where X± = Ix ± Hyi q = Kzz/e is the field gradient 
parameter, and rj = {Vxx — Vyy)/Vzz is the biaxial- 
ity (asymmetry) parameter. The former is the primary 
coupling constant of QI, and the latter determines the 
mixing between the free nuclear spin magnetic quantum 
numbers. In the most general case that we shall consider, 
in addition to quadrupolar part Hq, a nuclear spin will 
have interactions with dc and rf magnetic fields in the 
form "Hm = —jhl ■ Bq and ?/j.£ ~ —jhl ■ B'^^cosWj.ft, 
where 7 is the gyromagnetic ratio of the nucleus. The 
weak rf part can safely be treated perturbatively, whereas 
the dc magnetic field can be strong and gives rise to Zee- 
man effect; we shall abbreviate the stationary states un- 
der both quadrupolar and Zeeman splittings as the QZ 
spectrajS Individual matrix elements for obtaining the 
QZ spectra for an arbitrary EFG and Bq as well as the 
T/j-f-initiated transition rates are given in the Appendix 
section for the sake of completeness. 



is on) for which we use the value 0.035^^^ At a nu- 
clear site, R, V'(R-) is the electron wave function, that we 
simply approximate with a z- varying (height-dependent) 
Gaussian profile over the QD as V(p, z) cx e"^^''/^^^^^ , 
where D{z) is the z-dependent diameter of the truncated 
cone-shaped QD. 

The non-secular part of ?/jjf , which is responsible for 
the spin flip-flops as a direct process becomes energeti- 
cally impossible due to the presence of the external mag- 
netic field and the large energy mismatch between nu- 
clear and electronic Zeeman energies. On the other hand 
in the noncollinear case of major EFG principal axis be- 
ing tilted from the growth axis the secular part becomes 
much more interesting. When we express Iz operator in 
EFG-coordinate components we observe that this allows 
a nuclear spin transition without changing the spin ori- 
entation of the electron. There is still an energy cost for 
the nuclear spin flip but this is much less compared to 
non-secular term being on the order of a few neV, which 
is negligible compared to the spontaneous lifetime broad- 
ening of A; ~1 /ieV of the exciton state in the optically 
created electron spin pumping configuratioui^ There- 
fore, the rate of noncollinear secular hyperfine interaction 
(NCSHFI) can become significant. For a transition from 
state \i) to |j), this rate is given by 
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III. RESULTS 



A. Test structures 



(7) 



C. Noncollinear secular hyperfine interaction 

Quite commonly in the nuclear spin experiments there 
exists an electron in the QD with an optically oriented 
spin along the growth axisJ^"— Through the hyper- 
fine interaction it then polarizes the QD nuclear spins, 
termed as the dynamic nuclear polarization.^^ Deng and 
Hu have suggested that in the presence of quadrupolar 
mixing a new spin depolarization channel becomes possi- 
ble through the hyperfine coupling.— If we leave out the 
weaker anisotropic dipolar part in the s-type conduction 
band of III-V semiconductors,^^ the hyperfine interaction 
can be represented by the isotropic Fermi contact term^^ 




secular non-secular 

(6) 

where Sz , S± are the z-component and the rais- 
ing/lowering electron spin operators, Ajj£ is the hyperfine 
constant of the nucleus, is the average fraction of the 
time electron is inside the QD (i.e., hyperfine interaction 



There exists a plethora of different realizations for 
the compositional and structural profiles of the self- 
assembled QDs; for a very recent experimental review, 
see, Ref. |38|, and specifically for the InAs QDs, see, 
Ref. m. Guided by the samples used in recent nuclear 
spin experiments^^— we center our discussion around a 
QD of a truncated cone shaped with a base (top) diam- 
eter 25 nm (10 nm), placed over a 0.5 nm wetting layer, 
all fully embedded in a GaAs host lattice. The computa- 
tional supercell contains some 1,800,000 atoms, most of 
which for the host matrix region, with the InAs QD re- 
gion having 16,702 In and 15,432 As atoms. In the case of 
alloy mixing discussions, we randomly replace a fraction 
of the In atoms with Ga atoms. We start from a uni- 
form compressive strain in the QD region by setting its 
lattice constant to that of bulk GaAs. During the energy 
minimization we use periodic boundary conditions while 
allowing for the computation box to shrink and expand, 
so as to attain a zero pressure on the wallsj^ After this 
relaxation the InAs QD height settles to 2.93 nm. Our 
choice for such a high-aspect-ratio QD is again guided by 
the samples used in relevant experimentsji^r'^ In order 



4 



to extract the dependence on the aspect ratio and al- 
loy composition of the QDs, different heights and indium 
mole fractions are considered as well. 



B. Strain and quadrupolar splitting 

To have a broad overview, we would like to start with 
Fig. [T] displaying the atomistic profiles over (100) and 
(010) cross sections of InAs and Ino.7Gao.3As QDs. A 
compressive in-plane strain, e± = {cxx + (-yy)/'^ is seen to 
be mainly preserved after relaxation whereas for e^z it is 
released along the growth axis to its environs leaving the 
QD region with a tensile Czz value. The opposite signs 
for the in-plane and growth axis QD strain components 
is a reflection of the Poisson effect 1^ Thus, with the bi- 
axial strain defined as = — e_L, a compressive value 
follows directly from these two. It can be observed by 
comparing left and right panels in Fig. [1] that the com- 
positional variation docs not lead to a qualitative change 
on the strain behavior. In the case of quadrupolar energy 
parameter, vq the variation among different elements is 
substantial, while those within each element simply fol- 
lows that of the biaxial strain. 

To support these observations with more quantitative 
data, in Table I the strain and quadrupolar statistics of 
the InAs QD atoms are summarized. Since a common 
way of nuclear spin polarization is via that of an opti- 
cally oriented electron spin,— a relevant measure for the 
nuclear ensemble would be an electron wave function- 
weighted statistics. Hence, in Table I and other discus- 
sions we provide equal-weighted and electron envelope- 
square-weighted statistics for QD nuclear spins, taking 
the aforementioned Gaussian form for the latter. Start- 
ing with the strain components, In and As atoms look 
very similar. The exception to this is the shear strain 



InAs 



Ino.7Gao.3As 



measure that we define as £5 



which 



is significantly larger for As. This is caused by those As 
atoms on the interface forming heterobonds with In and 
Ga atoms. The shear strain is known to be crucial for the 
piezoelectric field, whereas in this context it also causes 
the tilting of the quadrupolar axes. Regarding the angu- 
lar behavior, first of all, the strain major principal axis 
lies almost on the in-plane direction, which is due to the 
dominance of e±. As a direct consequence of this strain 
profile, and in particular due to the relation Vzz = Sues, 
the EFG major principal axis is oriented very close to the 
growth direction. In contrast, for the V^x and Vyy EFG 
components, the strain components exx and eyy appear 
in opposite signs and largely cancel. A curious discrep- 
ancy between In and As occurs for the tilting angle of the 
major quadrupolar axis from the growth direction, Oz,q- 
This angle is about eight times larger for As than In. The 
tilting is caused by the shear EFG components, and the 
discrepancy between In and As mainly follows from the 
S44 coefficients, which are 10.0 x 10^^ statcoulomb/cm^ 
and 25.855 x 10^^ statcoulomb/cm'^, respectively."^^ In 
addition, as mentioned before, shear strain is also larger 
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FIG. 1. (Color online) Atomistic profiles®^ over (100) and 
(010) cross sections for InAs and Ino.7Gao.3As QDs. Bigger 
spheres correspond to In atoms. 



for As, and hence, the larger deviation of the As ma- 
jor quadrupolar axis from the growth axis. As an indi- 
rect consequence of this. As system has twice as large 
EFG biaxiality, i] — {Vxx — Vyy) /Vzz compared to 
In (cf. Table I). Another marked difference is about 
the quadrupolar energy gaps. If we temporarily as- 
sume a uniaxial case (77 — > 0), the energy difference be- 
tween m = ±1/2 and m — ±3/2 levels becomes hfQ = 
ie^qQI [2/(2/ — 1)]. This energy is more than twice as 
large for As (-8.3 MHz) compared to In (-4 MHz). The 
underlying reason is that In has / = 9/2 whereas As has 
3/2 nuclear spins. The Q values are however in favor of 
In, which has 0.860 x lO"^** cm^ compared to As having 
0.27 X lO"^"* cm^."^^ As a matter of fact, if one considers 
the full nuclear spin manifold energy span under QI, i.e., 
between m = ±1/2 and m = ±/, In has an extend of 
~36 MHz, which is about four times larger than that of 
As. If we now focus on the difference between envelope- 
square- and equal-weighted (in parantheses) results, the 
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TABLE I. The mean and standard deviation for certain atomistic quantities of the considered InAs QD. e^: in-plane strain 
perpendicular to growth axis, eg: shear strain, e^: biaxial strain, vq: quadrupolar (lowest) energy splitting, ry: EFG biaxiality 
parameter, 6^^^- polar angle the strain major principal axis makes with the growth axis, &z,q: polar angle the quadrupolar 
major principal axis makes with the growth axis. Envelope-square-weighted statistics are given, while the raw (equal-weighted) 
values are quoted in parentheses. 
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FIG. 2. (Color online) Monolayer-averaged axial variations of 
e_B Ezz, ^z,Q, and uq along the growth axis starting from the 
base of the QD without including the wetting layer. Filled 
(hollow) symbols are for InAs (Ino.yGao.a As) QD. The Gaus- 
sian envelope square weighting is used for the contribution of 
atomistic quantities within each monolayer. 



main deviations are seen to be on the shear strain, EFG 
tilt angles and biaxiality values, which indicate that there 
comes a substantial contribution from the inclined inter- 
face regions in the case of equal- weighted statistics. Also 
note that standard deviation for all quantities are higher 
when the interface region is included (i.e., under equal- 
weighted statistics). 

Next, in Fig. [2] we return to the effect of alloy mixing 
and examine the variation of monolayer-averaged strain 
and quadrupolar quantities along the growth axis, us- 
ing Gaussian envelope-square weighting. In all cases the 
quantities show a rather flat profile along the growth axis 
up until the top two monolayers. It can be observed that 




0.12 
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0.04 



Quantum Dot Aspect Ratio 

FIG. 3. (Color online) Aspect ratio (defined as the base diam- 
eter over the height) dependence of es and vq for an InAs QD. 
Base (top) diameter is kept at 25 nm (10 nm) while the QD 
height is varied. Envelope-square-weighted values are used. 



compared to pure InAs substitution of 30% of Ga 

changes the overall QD strain state in such a way that 
both e^z and are reduced, and in corollary, so is the 
energy splitting, i'q. On the other hand the tilt angle of 
the major quadrupolar axis from the growth axis, 9z,q is 
larger caused by the local random alloy configurations. 
A large 0z,q has a number of implications such as en- 
hancing the NCSHFI as will be discussed later. These 
observations are further supported quantitatively in Ta- 
ble II, which compares the statistics of e^, eg, es, vq, 77, 
0z.e, and Q for three different alloy compositions. With 
the reduction of indium composition from the pure InAs 
case, the in-plane and biaxial strain diminishes as well as 
the quadrupolar energy splitting, whereas, the increased 
random alloying enhances the EFG biaxiality and causes 
substantial deviation of the major strain and quadrupo- 
lar axes from the growth axis especially for the As and 
Ga nuclei. 

Finally, we consider the role of aspect ratio of the QD 
on the strain and quadrupolar state. To begin with, if 
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TABLE II. The effect of alloy composition on the Gaussian envelope-square- weighted statistics for the same atomistic quantities 
considered in Table I. Three different alloy compositions are compared for a In^Gai_a;As QD: x =1, 0.7, and 0.4. 







Indium 






Arsenic 




Gallium 




X =1 


X =0.7 


X =0.4 


X =1 


X =0.7 


X =0.4 


X =0.7 


X =0.4 


^± 


-0.057 


-0.046 


-0.024 


-0.056 


-0.045 


-0.018 


-0.044 


-0.015 


es 


0.004 


0.008 


0.008 


0.005 


0.011 


0.011 


0.008 


0.008 


es 


0.090 


0.067 


0.034 


0.086 


0.065 


0.031 


0.068 


0.033 


UQ (MHz) 


3.974 


3.846 


2.346 


9.804 


8.048 


1.677 


4.555 


2.037 


V 


0.042 


0.198 


0.431 


0.117 


0.414 


0.584 


0.239 


0.476 


0... 


89.3° 


87.7° 


82.3° 


88.0° 


86.3° 


74.5° 


88.3° 


70.5° 




1.7° 


3.2° 


12.4° 


12.9° 


17.9° 


41.5° 


5.3° 


17.4° 



we were to have a spherical QD embedded into a host 
of a different lattice constant, in the continuum approx- 
imation we would have only a hydrostatic strain and no 
biaxial strain and no quadrupolar shift. In accordance 
with this, as seen in Fig. [31 QDs with a large aspect ratio 
possess a larger biaxial strain causing large quadrupolar 
shifts. Note that due to the truncated cone shape, an 
anisotropy always remains regardless of the height of the 
QD. In self- assembled QDs, while the height is controlled 
to a very high precision, such as by capping and double 
capping techniques)^ this is not the case for the lateral 
dimension, which is essentially determined by the local 
growth kineticsJ^ Therefore, we can expect a significant 
variance on the quadrupolar shifts from dot to dot within 
the same sample, much like their light emission proper- 
ties. 
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FIG. 4. (Color online) QZ spectra of single In and As nuclei 
for the Faraday and Voigt geometries. 



C. NMR spectra 

Now, including a dc magnetic field, Bq, both Fara- 
day and Voigt geometries will be discussed where Bq is 
parallel and perpendicular to the (optical) growth axis, 
respectively. In Fig. |4] the QZ spectra of a single In and 
As nuclei are given. For the quadrupolar parameters we 
use those listed in Table I under the envelope-squared 
QD statistics. Some of the comparisons regarding the 
quadrupolar splittings between As and In were already 
made in the previous section. Here, we would like to focus 
on the evolution of the QZ spectra from the quadrupolar- 
to the Zeeman-dominated regime. One can easily realize 
the marked discrepancy between the Faraday and Voigt 
geometries: in the former, as the magnetic field is in- 
creased, the quadrupolar-split m — ^I state moves down 
in energy through a number of band crossings followed 
by a final anticrossing. The couplings are localized to 
the vicinity of these "crossings" . In the case of Voigt ge- 
ometry however, there are no such band crossings. The 
fundamental difference is illustrated in Fig. [S] Most im- 
portantly, QI being a rank-two tensor has a bilateral axis 
because of which there remains the ±to degeneracy in 
its spectrum, whereas the magnetic field is vectorial, i.e., 
unilateral. In the Faraday geometry, as the field increases 
the system loses its bilateral character while slightly ro- 



tating to align with the magnetic axis. In the Voigt ge- 
ometry, almost orthogonal setting of the two axes causes 
all the states to be mixed as the system evolves, as can 
be seen from the corresponding QZ spectra. 

We shall now make use of the QZ spectra in interpret- 
ing the NMR absorption spectra. To excite the nuclear 
spins a linearly polarized rf magnetic field perpendicu- 
lar to Bq is introduced. If EFG were both uniaxial (i.e., 
77 = 0) and its major principal axis were collinear with Bg 
then only a perpendicular rf magnetic field would cause 
transitions. In our QDs where these assumptions do not 
hold, we observed that an rf field parallel to Bq also cou- 
ples to nuclear spins yielding a similar spectra. However, 
their intensity is about two orders of magnitude weaker, 
and for this reason will not be included here. We assume 
that initially all nuclear spin states are equally populated 
giving rise to full spectra, which can be justified based 
on the much smaller energy separations among QZ states 
with respect to thermal energy at few Kelvins, as well as 
due to the presence of dipole-dipole interactions among 
the nucleii^i^ As a caveat, we would like to point out 
that some parts of the full spectra, to be discussed be- 
low, can become hindered depending on how the nuclear 
spin ensemble gets prepared under various experimental 
realizations. 
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Faraday Voigt 

FIG. 5. (Color online) Illustration of the magnetic field (Bo), 
and quadrupolar (Z) axes for Faraday and Voigt geometries. 
The single- and bi-cones represent the |/) and | ±/) states for 
the Zeeman and quadrupolar cases, respectively. 



The NMR spectra of InAs and Ino.7Gao.3As QDs are 
shown in Figs. [5] and [7] for Faraday and Voigt geometries, 
respectively. To assist their interpretation we also include 
those of the element-resolved and single nucleus spectra 
as well. As a matter of fact, the fingerprints of the single 
nucleus spectra can be readily identified on the overall 
QD cases. Not observing In nuclei in the Voigt geome- 
try NMR spectra, Flisinski et al. attributed this to the 
2.5 times smaller population of each spin state compared 
to a spin-3/2 nucleus.'''* However, our predictions show 
that In nuclei must have stronger rf absorption intensity 
compared to As, which stems from the 7/^ dependence 
in the rate expression (cf. Appendix); together with the 
aforementioned population discrepancy this results in an 
overall factor of about five in favor of In over As. For this 
reason. In nuclei are dominant on the total spectra. As 
marked in the Faraday geometry of single nucleus cases, a 
borderline between the quadrupolar and Zeeman regimes 
can be introduced, which corresponds to about LIT for 
the As and 1.5 T for the In nucleus. These can be taken 
as the effective magnetic fields, Bq up to which QI is 
strong. Even though it is not as distinct, one can ob- 
serve that same values also hold for the Voigt geometry. 
The transitions among m states are also labeled on the 
single As nucleus with its simpler spectrum, where we 
use on either side of Bq the asymptotic pure quadrupo- 
lar or pure Zeeman basis. For each case individually, the 
selection rule is^ Am = ±1, however, EFG biaxiality 
or the noncoUinear EFG axis with respect to Bo intro- 
duce higher-order transitions albeit with much weaker 
strength. 

The severity of the inhomogeneous broadening can be 
clearly observed from a comparison between the single 



nucleus and element-resolved spectra. As expected, the 
higher-order transitions, which are forbidden in the ab- 
sence of QI are highly broadened because of the atom- 
istic level strain variation over the QD to which they owe 
their existance. If we now focus on the overall spectra, we 
observe that for Faraday geometry, the so-called central 
transition, 1/2 ^ -1/2 is the sharpest among all, even in 
the presence of alloy mixing in luo.rGao.a As, with the rea- 
son being that for this transition QI has no influencei^ 
On the other hand for the Voigt geometry all of single- 
nucleus resonances are broadened due to orthogonality of 
the quadrupolar axis with respect to Bo. 

In Fig. [8] we compare the NMR spectra at 1 T and 5 T 
for the Faraday and Voigt geometries. Note that to gain 
a broader insight to the compositional effects we also in- 
clude the Ino.4Gao.6As case. Once again it can be verified 
that for the Faraday geometry the central transitions of 
all QDs coincide with the same sharp resonances, marked 
by three vertical arrows in the 5 T case, corresponding 
to As (36 MHz), In (47 MHz) and Ga (51 MHz) nu- 
clei. Another pivotal observation is that the Ino.7Gao.3As 
and Ino.4Gao.6As QDs have progressively narrower over- 
all spectral support with respect to InAs case. This effect 
of alloying reminds the motional narrowing in the case of 
atomic gases j^i^ however, in this case caused by quite a 
different reason, where the spin-9/2 manifolds of the In 
nuclei are partially replaced by the spin-3/2 manifolds of 
the Ga nuclei, with the latter having much narrower span. 
Therefrom this simply suggests that the NMR spectra, 
especially in the Faraday geometry can be utilized to ex- 
tract the indium mole fraction, which is one of the key 
unknown material parameters for a specific QD under 
consideration*^2i^ 



D. NCSHFI-mediated nuclear spin depolarization 

Finally, we discuss within the Faraday configuration 
the NCSHFI-mediated nuclear spin depolarization, or for 
that matter the polarization process as well, as it is also 
governed by the same matrix element (cf. Eq. [T]). We 
assume that nuclear spin is polarized through an elec- 
tron spin that is aligned along the growth axis (here, z- 
direction). The presence of the quadrupolar field with 
tilted principal axes introduces a complication, as to 
which QZ state is to be ascribed for such a polarized nu- 
clear spin. For this purpose, we introduce the so-called, 
maximally overlapping state, |imax}, selected from all 
QZ states i G {—/,...,/} that maximizes the overlap, 
(m^jz)!, where \mz) denotes a free nuclear spin state, 
which is aligned along the polarizing electron's spin di- 
rection. 

In the inset of Fig. [9] we see the evolution of the 
maximally- aligned state through a number of band 
"crossings" over the QZ states. The total out-transition 
time from such a state, hence the depolarization time, for 
an individual In or As nucleus under the mean quadrupo- 
lar field values (taken from Table I) are shown in Fig. [HI 
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FIG. 6. (Color online) The Faraday geometry NMR spectra for InAs and Ino.7Gao.3As QD nuclei (top row), together with their 
element-resolved contributions, As (center row), and In (bottom row) for the InAs QD (left panel, center and bottom rows), 
contrasted with respective single- nucleus spectra (right panel, center and bottom rows). 




FIG. 7. (Color online) Same as Fig[S]but for the Voigt geometry. 
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FIG. 8. (Color online) Comparison of NMR spectra at 1 T and 
5 T for InAs (dashed), Ino.7Gao.3As (solid), and Ino.4Gao.6As 
(dotted) QDs. Upper (lower) plot is for the Faraday (Voigt) 
geometry. The vertical arrows in the Faraday geometry for 
5 T mark the central transitions for In, As, and Ga nuclei. 



We observe that the depolarization is enhanced when the 
initial maximally- aligned state goes through band cross- 
ings with other QZ states, and the broader minimum 
occurs with the final band anticrossing. In the case of al- 
loy mixing (dashed lines) larger tilt of the EFG axis from 
the growth direction in general causes rapid depolariza- 




dc Magnetic Field (T) 



FIG. 9. (Color online) Nuclear (de)polarization time of indi- 
vidual In and As nuclei due to NCSHFI. QD mean EFG values 
are used with solid (dashed) lines for InAs (Ino.7Gao.3As). In- 
sets show the evolution of the maximally- aligned state (high- 
lighted) among all QZ states as a function of magnetic field. 
All plots are for the Faraday geometry. 



tion and the minimum magnetic field that this occurs also 
decreases. The minimum NCSHFI-mediated depolariza- 
tion is seen to be on the order of an hour.— Just as in 
the NMR spectra, above Bq level where Zeeman regime 
takes over, the NCSFHI gradually becomes weaker. 



IV. CONCLUSIONS 

The progress of the optically detected NMR techniques 
assure that the manipulation of a few-number-nuclei will 
not be a distant futurej^"— In line with this prospect, 
we present the strain and quadrupolar fields that these 
nuclei are exposed to at an atomistic level. First we 
summarize several structural and compositional under- 
pinnings: A high aspect ratio enhances the QI, and the 
interface regions introduce biaxiality and the tilting of 
the major quadrupolar principal axis from the growth 
axis. On the other hand, alloy mixing of gallium into the 
QD reduces both the strain and the quadrupolar energy 
splitting. The spectra for Faraday and Voigt geometries 
are quite distinct from each other. For the latter, all lines 
are inhomogeneously broadened due to the orthogonality 
of the quadrupolar axis with Bq. For the former, central 
transition, 1/2 — >■ -1/2 remains sharp even in the pres- 
ence of alloy mixing. Forbidden transitions are also ob- 
served, though highly broadened, arising from the EFG 
biaxiality and the tilting of the quadrupolar axis from the 
growth direction. The borderline between the quadrupo- 
lar and Zeeman regimes is extracted as 1.5 T for In and 
1.1 T for As nuclei. At this value the nuclear spin depo- 
larization rate due to the noncoUinear secular hyperfine 
interaction with a resident electron in the QD gets max- 
imized. In the case of alloy mixing larger tilting of the 
EFG axis from the growth direction causes more rapid 
depolarization and the minimum magnetic field that this 
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occurs also decreases. The shortest NCSHFI-mediated 
depolarization time is seen to be on the order of an hour. 
As Zeeman regime takes over above Bq, this depolariza- 
tion channel progressively becomes weaker. 

Note added: After the initial submission of this work 
a new ODNMR technique is reported?^ the experimen- 
tal data for InGaAs QDs at 5.3 T agree extremely well 
with the central transitions marked in Fig. El with the 
added feature that the less common ^^Ga isotope is also 
resolved, whereas in this work we consider only the dom- 
inant ^^Ga isotope. 



of the parts of the full Hamiltonian as 



{m\HM\m) 



1 



mcosc 



-I- - (sin 6* cos ( 



{m'\nQ\m) = Aq { [3m2 - /(/ + 1)] S^'^m 

Tj ^ 



where fi{m) = //(— m — 1) = y^{I — ■m){I + m + 1), and 
= 41(21% ^ Solving for the QZ spectrum 

essentially yields the expansion coefficients, of the 
QZ states \i) in terms of free spin states as 
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APPENDIX: MATRIX ELEMENTS AND RF 
TRANSITION RATES 



iz) = 



(8) 



771— — I 



In the case of an incident rf field, the nuclear spins are 
excited over their established QZ spectrum through the 
Hamiltonian 



n 



rf - 



Bx^x + BylY + B^^Xz 



cos Wj-fi 



(9) 



Hence, based on Fermi's golden rule the rf absorption 
rate from an initial state \i) to any final state |j) will be 



2A/h 



(10) 



where A is the fundamental linewidth of an individual 
nuclear spin for which we take 10 kHz for all the nucleus 
types in this work.^*- The corresponding matrix element 
is then given by 



In the orthogonal coordinate system as defined by the 
principal axes of the EFG, the dc magnetic field vector 
Bo will be in general tilted as governed by the spherical 
polar angles 6, and (j) so that its contribution becomes^ 



n 



M 



-Ml {Ix sin 6 cos + ly sin sin cf) + Iz cos ( 



where — jBq. Choosing the angular momentum quan- 
tization axis as Z, and denoting^'^ the free nuclear spin 
states by |to), we can easily obtain the matrix elements 



m—^I 



B^^ (C^^_i) Cl^Ji{- 



where B^f = (s^f ± B'f^ /2. In the context of NCSHFI 
having essentially the same matrix element one 
has to replace the components of B'^'' in the above ex- 
pression with the components of the unit vector along 
the growth axis, z expressed in the EFG coordinate sys- 
tem. 
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